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ABSTRACT 

We present initial results from a new high-contrast imaging program dedicated to stars that 
exhibit long-term Doppler radial velocity accelerations (or "trends" ) . The goal of the TRENDS 
(TaRgetting bENchmark-objects with Doppler Spectroscopy and) imaging survey is to directly 
detect and study the companions responsible for accelerating their host star. In this first pa- 
per of the series, we report the discovery of low-mass stellar companions orbiting HD 53665, 
HD 68017, and HD 71881 using NIRC2 adaptive optics (AO) observations at Keck. Follow-up 
imaging demonstrates association through common proper-motion. These co-moving compan- 
ions have red colors with estimated spectral-types of K7-M0, M5, and M3-M4 respectively. We 
determine a firm lower-limit to their mass from Doppler and astrometric measurements. In the 
near future, it will be possible to construct three-dimensional orbits and calculate the dynamical 
mass of HD 68017 B and possibly HD 71881 B. We already detect astrometric orbital motion 
of HD 68017 B, which has a projected separation of 13.0 AU. Each companion is amenable to 
AO-assisted direct spectroscopy. Further, each companion orbits a solar-type star, making it 
possible to infer metallicity and age from the primary. Such benchmark objects are essential for 
testing theoretical models of cool dwarf atmospheres. 

Subject headings: keywords: techniques; radial velocities, image processing, high angular resolution; 
astrometry; stars: individual (HD 53665, HD 68017, HD 71881), binaries, low-mass, brown dwarfs 
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1. INTRODUCTION 

It is commonly thought that the radial ve- 
locity (RV) method for detecting companions to 
nearby stars only provides information about low- 
mass bodies with short orbital periods. However, 
Doppler measurements are quite sensitive to dis- 
tant objects, because the RV semi-amplitude, K, 
decreases slowly with orbital period, P (or semi- 
major axis, a), according to: 



K oc m sin(i) P"^/^ 
(X TO sin(i) a^^^'^ 



(1) 
(2) 



where to is the companion mass and i is the or- 
bit inclination. With ^1-10 m/s precision, M- 
dwarfs, brown dwarfs, and super- Jupiters are de- 
tectable out to tens of AU, it just takes a long 
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time t o complete a full orbit (e.g., [Howard et al 
Despite having knowledge of only a frac- 
tion of an orbital cycle, RV accelerations (trends) 
are tremendously useful: they show us conclu- 
sively that something initially hidden from view 
is tugging on the visible star. 

Informed target selection is an order of mag- 
nitude effect for substellar companions. For 
example, wide-separation (a > 10 AU) brown 
dwarfs are rare, having been found to orbit only 
sa3% of solar- type stars ( Metchev and Hillenbrand 
20091 ). And, while companions in the planetary- 
mass regime are expected to be more common, 
an upper-limit to the frequency of super- Jovian 
(m > 3Mj) bodies over a similar semi-major 
axis range is set at ~20% for solar-type stars 
(iNielsen and Closdl2010l: Iviean et al.ll2012l) . By 



observing a sample of intrinsically companion- 
rich sources, those with clear Doppler accelera- 
tions, it is possible to by-pass the inefficiencies 
common to high-contrast programs that nomi- 
nally select sta rs based solely on age and proxim- 



2007: Lafreniere et al. 


2007; .Leconte et al. 


.2010f 


Ehrenreich et al. 


20101). 



In addition to high observing efficiency, there 
are significant scientific benefits to combining 
Doppler measurements with high-contrast imag- 
ing. RV and direct astrometry observationtlj 
may be used in concert to calculate the compan- 
i on orbit (all 6 el ements) and dyn amical mass 
(iBoden et al.l bood) Fl For example, ICrepp etlir 



20121 have measured the mass of the benchmark 



brown dwarf HR 7672 B with a fractional uncer- 
tainty of 4% by monitoring its motion over ^33% 
of an orbit cycle. Masses determined independent 
of photometry and spectroscopy inform theoret- 
ical atmospheric models by helping to break de- 
generacies between the various input parameters, 
such as mass, radius, age, e ffective temperature, 



and chemical com position ([Barman et al.l 12011 



Janson et al.l 120111) . Dynamical masses are like- 



^ By direct astrometry we are referring to following the com- 
panion along its orbit with direct imaging, measuring the 
sky-projected separation and position angle relative to the 
star as a function of time. Indirect astrometry involves 
measuring the position of the star relative to other (dis- 
tant) nearby stars or fiducial reference points. 

^See lRodigas et al.ll20lil for a discussion regarding mass and 
orbit constraints for the case of Doppler accelerations com- 
bined with direct imaging non-detections. 



wise important for calibrating thermal evolution- 
ary models, providing a measure of subste llar ob- 
jects luminosity as th ey fad e with time (St evenson 
19911: iBurrows et al]|l997t ). Further, if the com- 



panion orbits a solar-type star, the mctallici ty and 



ICrepp et alboi^ lBiller et al.ll2010h ." 



age may be inferred from the prima ry (L iu et al 
' 20071: iDupuy et al.|l2009l: Ijohnson and App^ liooa 



With this motivation in mind, we have estab- 
lished an interdisciplinary program that uses a 
combination of Doppler observations and high- 
contrast imaging. Specifically, we use years of pre- 
cise RV measurements to identify promising tar- 
gets for follow-up high-contrast observations. The 
goals of our program are to: 

1. detect companions with mass in the m « 
5M,j - 500M,7 range; 

2. acquire spectro-photometric measurements 
across a wide bandpass (YJHKLM); 

3. perform follow-up direct astrometric mea- 
surements to break the sin i inclination de- 
generacy resulting from Doppler measure- 
ments and calculate dynamical masses; 

4. test theoretical atmospheric and evolution- 
ary models in a regime where they currently 
break-down (low temperatures). 

Based primarily at Keck Observatory, our RV 
measurements are obtained using the Hig h Reso 



lution Echelle Spectrometer (HIRES; iVogt et al 



19941 ) at Keck I. These measurements are of- 
ten augmented by previous and concurrent ob- 
servations at Lick Observatory. High-contrast 
imaging observations are obtained using NIRC2 
(PI: Keith Matthews) and th e Keck II AO sys- 



tem (jWizinowich et al.l l2000| ). We are also ex- 



panding this program to include AO observations 
with the Project 1640 s pectral-imager at Palo- 
mar (jHinklev et al.ll201l[ ) and LMIRCam at the 
Large Binoctilar Te lescope ( Skrutskie et al. 2010t 



Skemer et al.l l2012f) . Our strategy is to maxi- 
mize on-sky sensitivity by employing all of the 
powerful techniques recently developed for high- 
contrast imaging applic ations, including coron- 
agraphy (Guvon et ah ,200l ICrepp et al.] l2010l ) 
and aggressive point-spread function subtraction 
to re move residu al scattered starlight from im- 
ages (Marois et al. 2006; iLafreniere et al.i i2007i 
Crepp et al.ll201ll ). 
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We have found that selecting targets based on 
the presence of an RV trend is an effective ap- 
proach, allowing one to take a "short-cut" for 
finding stars likely to host a low-mass compan- 
ion amenable to direct imaging detection. In this 
paper, we report the the discovery of three M- 
dwarf companions orbiting solar-type stars. Each 
companion has red colors and a low luminosity, 
otherwise they would be noticed as double-lined 
spectroscopic binaries in RV data. The compan- 
ions are amenable to AO-assisted sp ectroscopy 
(|Bowler et al.ll2010l : IPuevo et aLll2012[ ) and repre- 
sent the first discoveries of our program. 

2. OBSERVATIONS 

2.1. High-Resolution Stellar Spectroscopy 

Doppler Measurements 

Precise RV data were obtained with t he High- 
Reso lution Echelle Spectrometer (HIRES jVogt et al 
19941 ) at Keck. We use the iodine cell refer- 
encing method to calibr ate instrument drift and 
measure Doppler shifts (jMarcv and Butleiill992 : 



Butler et al. 19961 ). Observations for HD 53665 



HD 68017 and HD 71881 began on Jan. 25, 1998, 
Jan. 13, 1997, and Dec. 24, 1997 respectively In 
each long-term acceleration indicated that 

the star is orbited by a distant body (Fig. 1). Ac- 
celerations are approximately linear for HD 53665 
and HD 71881, whereas HD 68017 shows signifi- 
cant orbit curvature (change in the acceleration). 
In ^4.21 we use Doppler measurements in combi- 
nation with imaging observations to constrain the 
mass of each companion. 

Stellar Properties 

Stellar (template) spectra, taken with the io- 
dine gas cell removed from the optical path, 
were analyzed using the LTE spectral synthesis 
code Spe ctroscopy Made E asy (SME) described in 
Valenti an d 

Fische] l2005i SME provides an es- 
timate of the stellar effective temperature (Tcff), 
surface gravity (logg), metaUicity ([Fe/H]), and 
projected rotational velocity (wsini). Table [5] 
lists the spectral-type and physical properties of 
each star derived from spectral fitting along with 
comparison to theoretical isochrones. 
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Fig. 1. — Precise Doppler radial velocity measure- 
ments. Our time baseline exceeds a decade for 
each star. HD 68017 shows significant curvature. 
It will be possible to calculate the dynamical mass 
of HD 68017 B with only several more astrometric 
measurements. 



2.2. Adaptive Optics Imaging 

HD 53665, HD 68017, HD 71881 were each dis- 
covered on the same night, February 22, 2011 UT, 
using NIRC2 at Keck. We used the narrow cam- 
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Fig. 2.— Discovery images of HD 53665 B, HD 68017 B, and HD 71881 B taken in the K' filter (A^ = 2.12 ^m) 
with NIRC2 on Feb. 22, 2011. Companions are indicated by an arrow. 



era setting (9.963 ± 0.006 mas pix~^ plate-scale 
(jGhez et al.l l2008f ) ) to provide fine spatial sam- 
pling of the system point-spread function. Images 
were initially obtained in the K' (Ac — 2.12/im) 
filter for search-mode operation. Companions to 
each star were noticed in raw frames. We executed 
a three-point dither pattern to facilitate removal 
of sky background noise. The companion orbiting 
HD 68017 is fainter and closer to its host star com- 
pared to the other targets, so we also obtained im- 
ages placing the star behind the (partially trans- 
missive) 300 mas diameter coronagraphic mask. 

Images were processed by flat-fielding, correct- 
ing for hot pixels with interpolation, subtracting 
the sky background, and rotating the frames to 
standard north-east orientation. Fig. 2 shows fully 
processed images of each companion. We acquired 
follow-up observations on January 7, 2012 UT in 
different filters to obtain colors and assess whether 
the candidates were associated with their host 
star. The angular separation and position angle 
of each companion are listed in Table 1. Bright- 
ness ratios are listed in Table 2. 

3. ASTROMETRY 

Our astrometric observations consist of two 
epochs separated by 0.9 years for each source 
(Table 2). All three stars have large proper- 
motions (see Table 1), allowing us to easily de- 
termine whether the companions share the same 
space motion over this time frame. We measured 



an accurate separation and position angle of each 
companion follow ing the technique described in 
Crepp et al.l l2012i We first fit two-dimensional 



Gaussian functions to the stellar and companion 
point-spread functions to locate their centroids in 
each frame. The primary star was not saturated 
in any of our dithered images. We then correct 
for distortion in the NIRC2 focal plane using pub- 
licly available solutions provided by Keck Observa- 
tory's astrometry support paged. The results are 
averaged and uncertainty in the separation and 
position angle is taken as the standard deviation, 
taking into account uncertainty in the plate scale 
and orientation of the array by propagating these 
errors to the final calculated position. 

Fig. 3 shows multi-epoch astrometry measure- 
ments plotted against the expected motion of 
a distant background object. We find that all 
three companions, HD 53665 B, HD 68017 B, and 
HD 71881 B, are each clearly associated with their 
respective host stars, implying that they are grav- 
itationally bound. The positions of HD 53665 B 
and HD 71881 B have changed by R:!l-2 pixels 
over the 0.9 year time-frame. Meanwhile, an un- 
related background source placed at infinite dis- 
tance would have moved relative to the host star 
by 163.7 mas (16.4 pix) and 296.5 mas (29.8 pix) 
respectively. HD 68017 B has the smallest pro- 
jected separation of the three (13.0to;2 ^U, Feb. 
2011) and appears to exhibit significant orbital 



^ http: / / www2 . keck, hawaii . odu /inst / nirc2 / forReDoc / post_observing/dewarp / 
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Target Name Date (UT) JD-2,450,000 p (mas) P. A. (degrees) Proj. Sep. (AU) 



HD 53665 B 



HD 68017 B 



HD 71881 B 



Feb. 22, 2011 
Jan. 7, 2012 
Feb. 22, 2011 
Jan. 7, 2012 
Feb. 22, 2011 
Jan. 7, 2012 



5614.79 
5933.93 
5614.81 
5933.96 
5614.81 
5933.97 



1420.7 ± 1.0 
1413.4 ±0.9 

594.5 ±0.5 

574.6 ±0.5 
851. 7± 1.1 
859.3 ±0.6 



21.3° ±0.1° 
20.9° ±0.1° 
248.2° ±0.1° 
240.3° ±0.1° 
247.8° ±0.1° 
246.7° ±0.1° 



102.9 
102.4 



6.0 
+5.0 
6.1 

c+O.l 
3-0.2 
2+0.9 



12. 
35. 
35. 



1.1 

C-+0.9 
0-1.1 



Table 1: Summary of astrometrie measurements. 
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Fig. 3. — Astrometrie measurements (red-crosses) demonstrating that each companion is comoving with its 
host star. Sohd-curves show the path that a distant background object would execute over the same time 
frame accounting for proper-motion and parallactic motion from Feb. 22, 2011 through Jan. 7, 2012. 



motion in a clockwise direction. 

4. COMPANION MASS ESTIMATE 
4.1. Mass from Photometry 

We measured the brightness of each com- 
panion relative to its host star by performing 
aperture photometry, accounting for contami- 
nation from the primary. Stellar magnitudes 
were first converted fro m 2MASS measurements 
(ISkrutskie et all l2006bh to the MKO filter sys- 
tem (^Tokunaga et al. 'j2002h using transformations 
from Carpenter 200 lO Our observations were ac- 
quired in the J,H,K' (MKO) bands. We have 
assumed that Aif k, Aif', an assumption that is 
justified given the relatively mild colors of solar- 
type stars, and fact that uncertainty in the mea- 
sured magnitude difference and parallax domi- 
nate the uncertainty in absolute magnitude. The 



mass of each companion is found by comparing its 



brightness to late-type dwarfs using iDotter et al 



[2OO8. theoretical evolutionar y tracks (the Dart- 
mouth models), and also the Delfosse et al. 2000l 
empirical relations which correlate Mk with mass. 
Differential, apparent, and absolute magnitudes 
are listed in Table 2. 

We find that HD 53665 B has a mass of 0.65 ± 
O.O3M0 based on photometry (jDotter et al.ll2008f) . 
Its absolute magnitude in each near-infrared band 
is consistent with this value . Using Table 5 from 



Kraus and Hillenbrand 2007. HD 53665 B has col 



ors and brightness consistent with either an K7 
or MO dwarf. For comparison, we find a mass 
of 0.70 ± 0.03 Mf:^ based on the measured Mk' 



4.81 ± 0.16 (jDelfosse et al.ll2000^ 



*http:/ /www. astro. caltech.edu/ jmc/2mass/v3/transformations/ 



HD 68017 B has a mass of ().16± 0.02M(t, based 
on photometry ( Dotter et al.l 12008 ). Its absolute 
magnitude is estimated only in the K-band, as 
our unocculted H-band observations of the pri- 
mary star from Jan. 7, 2012 were saturated. Ta- 
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HD 53665 A 


HD 68017 A 


HD 71881 A 


R.A. (J2000) 


07 05 52.8 


08 11 38.6 


08 31 55.0 


Decl. (J2000) 


-01 01 13.7 


+32 27 25.7 


+50 37 00.1 


B (mag) 


7.76 ±0.02 


7.50 


8.06 


V (mag) 


7.26 ±0.01 


6.81 


7.43 


^2MASS (mag) 


6.294 ±0.024 


5.48 


6.284 + 0.020 


-ff2MASs (mag) 


6.066 ±0.027 


5.15 


6.027 ±0.027 


-f'^2MASS (mag) 


5.983 ±0.018 


5.09 


5.959 ±0.029 


d (pc) 


79 9+3.9 
'^•^_4.0 


21.8^°- 


41.3^}^ 


p.m. (mas/yr) 


9.4, -15.1 


-462.6, -644.2 


-81.4, -338.6 


Spec. Type 


F8V 


G4V 


GIV 


M,(Mq) 


1 CI +0.20 
-•^■^^--0.27 


85+0-04 

^■°3-0.03 


1 04+0-06 


Tcff' (K) 


6225 ± 44 


5552 ± 44 


5821 ± 44 


log g (cm/s^) 


4.05 ±0.06 


4.65 ±0.06 


4.29 ±0.06 


[Fe/H] 


0.17 ±0.04 


-0.44 ±0.03 


-0.05 ±0.03 


vsini (km/s) 


8.6 ±0.5 


0.8 ±0.5 


2.1 ±0.5 




-5.022 ±0.009 


-4.928 ±0.013 


-5.043 ±0.002 


Age (Gyr) 


9 1 +0.8 
^■-'^-0.4 




4 3+1-0 





HD 53665 B 


HD 68017 B 


HD 71881 B 


AJmko 


3.85 ±0.15 






Ai/MKO 


3.26 ±0.07 


> 4.16 


4.26 ±0.04 




3.14 ±0.10 


4.92 ±0.10 


4.12 ±0.08 


Jmko 


10.12 ±0.16 






Hmko 


9.31 ±0.08 


> 9.29 


10.26 ±0.07 


Kmko 


9.11±0.11 


10.00 ±0.11 


10.07 ±0.09 


-^Jmko 


5.82 ±0.20 






^Hmko 


5.01 ±0.14 


> 7.60 


7.18 ±0.09 


^K'mko 


4.81 ±0.16 


8.31 ±0.11 


6.99 ±0.11 


Spec. Type 


K7-M0 


M5 


M3-M4 


rridyn {Mq) 


> 0.63 


> 0.08 


> 0.17 


TOmodcl (Mq) 


0.65 ±0.03 


0.16 ±0.02 


0.31 ±0.03 


'Tlempirical (Mq) 


0.70 ±0.03 


0.15 ±0.01 


0.29 ±0.02 



Table 2: (Top) Coordinates, apparent magnitude, distance, proper motion (p.m.), spectral-ty pe, and physical 
proper ties of each host star. Near-infrared magnitudes for the primary star are from 2M ASS ( Skrutskie et alj 



2006af ). Distance estimates are based on measured parallax from the Hipparcos satellite (jvan L ccuwcn 2007|). 
Effective temperature (Tes), surface gravity (log g), and metallicity ([Fe/H ]) are derived using SME Ages 
are estimated using theoretical isochrones for HD 53665 {B — V ~ 0.47 ) (IValenti a nd Fischei"2005'), and 
gyrochronology for HD 68017 (B-V = 0.69) and HD 71881 {B-V = 0.63) (jMamaiek and Hillenbrand 200d). 
(Bottom) Companion magnitude difference, absolute magnitude, est imated spectral-ty pe, mass constraint 
from dynami cs (mdim) and estima ted mass from photometry using the lDotter et al.]|2008l atmospheric models 
(wmodei) and Delfosse et al. 2000l M^-mass empirical relations (mempiricai)- 



ble 2 shows a lower-limit for the H-band magni- 
tude. HD 68017 B has a K-band brightnes s con- 
sisten t with an M5-dwarf ( Kraus and Hillenbrand! 
20071 ). For comparison, we find a mass of 0.15 ± 
O.OIMq based on the measured Mk' = 8.31 ±0.11 



(jPelfosse et al.ll2000f) . 

HD 71881 B has a mass of 0.31 ± O.O3M0 
based on photometry. Its absolute magnitude in 
the H and K bands are both consistent with this 
value. HD 71881 B has colors and brightness 
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most consistent with either a n M3 or M4 dwarf 
( Kraus and Hillenbrand! 2007 ). For comparison, 
we find a mass of 0.29ib0. 02Mf:^ based on the m ea- 
sured Mk' = 6.99 ±0.11 (|Delfosse et al.ll2000[) . 



4.2. Mass Lower-limit from Dynamics 

When combined with a RV trend, a single epoch 
of imaging observations (single measurement of 
the projected physical separation) places a lower- 
limit on the compan ion dynamical mass ( Torres! 
19991 : iLiu et al.l[2002h . We have measured the in- 
stantaneous Doppler acceleration for HD 53665 
and HD 71881, assuming the RV data may be ap- 
proximated as linear across the full time baseline. 
Using a Markov chain Monte Carlo (MCMC) anal- 
ysis, we find slopes of: 



('^«/'^0hD53665 = 
('^«/'^0hD71881 = 



+5.3±0.3ms"Vr"^ 
-10.3 ± 0.2 ms" Vr"^- 



With a projected separation of 102.9 AU (Feb. 
22, 2011), HD 53665 B has a minimum dynam- 
ical mass of 0.77 ± O.14M0. Accounting for un- 
certainty in the stellar parallax, measured angular 
separation, and RV acceleration, we adopt a min- 
imum mass of O.63M0. This value is consistent 
with the mass derived from photometry provided 
HD 53665 B has an edge-on orbit (modulo faulty 
assumptions regarding system coevolution, or sys- 
tematic errors in the isochrone models). Likewise, 
HD 71881 B has a projected separation of 35.2 
AU (Feb. 22, 2011), which corresponds to a min- 
imum dynamical mass of O.17M0. Constraints 
from Doppler RV and imaging data are compared 
to mass estimates from photometry in Table 2. 

HD 68017 B has a projected separation of 
13.0tg;2 AU (Feb. 22, 2011). Evaluating the 
local RV slope using the most recently obtained 
21 data points (from Sept. 25, 2010 to present), 
we find (dwM)„„^3„^, = +16.3 ± 0.9 ms-iyr-\ 
which corresponds to a minimum dynamical mass 
of O.O3M0. However, the full RV time series 
shows significant curvature (change in the accel- 
eration), enabling a more sophisticated analysis. 
We have performed MCMC simulations that si- 
multaneously fit the Doppler and astrometric data 
using a Keplerian orbit. Our calculations account 
for the measured orbital motion of the companion 
from both imaging epochs. While the currently 
available data set provides insufficient information 



to calculate a unique orbit inclination, we find that 
the companion minimum mass constraint becomes 
mdyn > O.O8Af0 (68.2% confidence). 

5. SUMMARY 

We have established a new observing program 
that uses precise RV measurements to identify 
promising targets for high-contrast imaging obser- 
vations. The nearby, solar-type stars HD 53665, 
HD 68017, and HD 71881 exhibit long-term 
Doppler accelerations ("trends"). We have used 
NIRC2 AO observations at Keck to directly de- 
tect the companions responsible for causing the 
trend in each case. Follow-up astromctry demon- 
strates that each candidate is co-moving with the 
primary star. Relative photometry measurements 
suggest spectral-types of K7-M0, M5, and M3-M4 
respectively. 

As inferred from their host star, each M- 
dwarf companion is a metallicity and age bench- 
mark object. With continued Doppler monitor- 
ing and follow-up AO observations, two out of 
the three companions, HD 68017 B and possi- 
bly HD 71881 B, will also serve as mass bench- 
marks. HD 68017 has a projected separation of 
only 13.01q2 and already shows significant RV 
curvature (change in the acceleration) and mea- 
surable astrometric motion, making it possible to 
estimate a dynamical mass in the next several 
yearsH Using currently available data, we calcu- 
late a lower-limit to the mass of each companion. 
We posit that HD 53665 B must have a near edge- 
on orbit given the agreement between the mass 
estimate from photometry and dynamics. 

Few mass, age, and metallicity benchmark 
dwar fs are known to date (iLiu et al. l2007t iDupuv et al 



20101: iBowler et all |2012a||b|). The goal of the 



TRENDS high-contrast imaging program is to 
discover and characterize low-mass stellar and 
substellar companions with physical properties de- 
termined independently from spectro-photometric 
measurements, in order to calibrate theoretical 
atmospheric models and thermal evolutionary 
models. Each of the companions presented are 



^For comparison, given the typical (high) signal-to-noise ra- 
tio of Doppler measurements for solar- type stars, and direct 
astrometric measurements, iCrepp et al.ii2012l have shown 
it possible to calculate accurate dynamical masses for com- 
panions with semimajor axes as large as Ril9 AU. 
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amenable to direct spectroscopy using AO-fed 
integral-field units, such as OSIRIS at Keck or 
Project 1640 at Palomar. 
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Table 3 

Radial Velocities for HD 53665 



HJD 


RV 


Uncertainty 


-2,440,000 


(m s-i) 


(m s-i) 


10838.8830 


-25.83 


5.69 


10861.8568 


-9.69 


5.98 


11071.1199 


-16.24 


5.84 


11171.9255 


-14.63 


5.79 


11226.8781 


-17.12 


5.65 


11544.0550 


-1.81 


5.92 


11552.9643 


-20.22 


6.06 


11585.9500 


-23.16 


5.82 


11883.0298 


-33.55 


5.82 


11973.8361 


-22.17 


5.91 


12235.9597 


-16.01 


5.60 


12574.0727 


0.74 


5.70 


12600.9894 


-7.03 


5.60 


12653.0320 


15.80 


5.83 


12987.9295 


-13.35 


5.77 


13304.1114 


16.06 


5.24 


13400.9519 


22.45 


5.38 


14024.1377 


31.39 


5.45 


15521.9975 


32.68 


5.38 


15672.8256 


57.91 


5.50 


15871.0233 


43.77 


5.25 
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Table 4 

Radial Velocitier for HD 68017 



HJD RV Uncertainty 



-2,440,000 








10461 


, u 


50.69 


2, 


.36 


1 0784 




41.84 


2, 


.41 


1 0807 




30.26 


2, 


.74 


1 08*^7 




39.68 


2, 


.46 


1 08*^8 




35.12 


2, 


.53 


1 0861 


.8262 


39.34 


2, 


.46 


1 0869 


7^86 


42.87 


2, 


.38 


11171, 




18.20 


2, 


.44 


11226, 


7940 


31.60 


2, 


.51 


11227, 




30.02 


2, 


.43 


11229, 


.9098 


32.27 


2, 


.33 


11551, 


.0555 


17.13 


2 


.61 




91 88 

. t7 -LOO 


15.21 


2, 


.46 


11884, 


.1426 


14.94 


2 


.52 


1 1 8Q8 


.1321 


21.48 


2 


.50 


1 1 8QQ 


1 905 


23.57 


2, 


.41 


1 1 QOO 


0979 


18.66 


2 


.33 


1 1 Q01 


1 5*^7 

. J- (JO i 


16.07 


2 


.40 


11972, 


.0135 


11.99 


2, 


.44 


1 1 Q72 


9874 


15.01 


2, 


.39 


1 1 97*^ 


8854 


16.30 


2 


.45 


1 1 974 


8699 


12.93 


2 


.49 


12243, 


0488 


5.48 


2, 


.58 


iZO / O, 


.iUUD 


2.04 


2 


.53 


12601, 


.0207 


-2.07 


2, 


.44 


12680, 


.9684 


-8.36 


2 


.47 


12711, 


.7407 


-2.32 


2 


.58 


12988, 


.0081 


-12.00 


2, 


.55 


13303, 


.1395 


-6.03 


2 


.78 


13398, 


.8550 


-15.14 


2 


.12 


13692, 


.9970 


-14.80 


2, 


.13 


13693, 


.0957 


-16.00 


2, 


.12 


13694, 


.1285 


-17.28 


2 


.13 


13695, 


.1055 


-16.35 


2, 


.13 


13696, 


.0847 


-16.56 


2 


.19 


13697, 


.1010 


-28.71 


2 


.18 


13724, 


.0212 


-17.79 


2 


.15 


13725, 


.0429 


-18.73 


2 


.15 


13747, 


.0478 


-12.93 


2 


.14 


13747, 


.9616 


-15.76 


2 


.15 
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Table 4 — Continued 



HJD RV Uncertainty 



-2,440,000 


(m s "'") 




(m s "'") 


13748, 


.9262 


-16.15 


2, 


.37 


13749, 


.8534 


-13.49 


2 


.16 


13750, 


.8602 


-15.83 


2, 


.15 


13751, 


.9266 


-10.32 


2 


.39 


13752, 


.9635 


-14.96 


2 


.14 


13753, 


.9402 


-12.69 


2, 


.13 


13775, 


.7943 


-13.93 


2, 


.14 


13776, 


.9239 


-10.81 


2 


.14 


14129, 


.9829 


-16.35 


2, 


.14 


14399, 


.1317 


-10.94 


3 


.11 


14492, 


.9688 


-22.14 


2 


.43 


14547, 


.9008 


-16.94 


2, 


.46 


14548, 


.8355 


-11.53 


2 


.53 


14809, 


.9817 


-6.87 


2 


.59 


14867, 


.8982 


-8.70 


2, 


.48 


14927, 


.8853 


-12.18 


2, 


.56 


14929, 


.8538 


-2.49 


2 


.48 


14984, 


.7583 


-8.21 


2, 


.42 


15109, 


.1438 


-5.35 


3 


.02 


15134, 


.1513 


0.84 


3 


.02 


15164, 


.0166 


1.78 


2, 


.44 


15172, 


.1149 


-4.52 


2, 


.52 


15188, 


.0140 


-2.65 


2 


.51 


15189, 


.9820 


-1.62 


2, 


.44 


15192, 


.0018 


-0.65 


2, 


.49 


15198, 


.9722 


-8.83 


2 


.50 


15229, 


.1011 


-1.10 


2, 


.16 


15229, 


.7942 


-5.34 


2, 


.17 


15231, 


.8098 


2.19 


2 


.17 


15251, 


.9199 


0.16 


2, 


.17 


15255, 


.7629 


-5.62 


2 


.15 


15260, 


.7790 


-6.72 


2 


.16 


15289, 


.7242 


2.54 


2, 


.50 


15311, 


.7447 


3.50 


2 


.49 


15344, 


.7419 


-6.89 


2 


.69 


15465, 


.1494 


5.98 


2, 


.36 


15468, 


.1055 


4.81 


2 


.36 


15469, 


.1528 


3.65 


2 


.37 


15470, 


.1502 


9.43 


2 


.36 


15487, 


.1576 


5.54 


2 


.35 
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Table 4 — Continued 



HJD RV Uncertainty 

-2,440,000 (m s"!) (m s"!) 



15490, 


.1555 


8.16 


2 


.41 


15491, 


.1524 


8.56 


2 


.35 


15501, 


.1589 


9.99 


2 


.38 


15522, 


.0236 


2.12 


2, 


.27 


15543, 


.1041 


3.73 


2 


.16 


15605, 


.9635 


16.44 


2 


.19 


15634, 


.8186 


8.89 


2, 


.16 


15671, 


.7972 


14.31 


2, 


.13 


15697, 


.7307 


12.49 


2 


.17 


15698, 


.7813 


11.57 


2, 


.18 


15843, 


.0775 


24.08 


2 


.38 


15880, 


.1365 


28.22 


2 


.93 


15904, 


.1782 


22.53 


2, 


.58 


15960, 


.9100 


33.29 


2, 


.55 


15999, 


.7837 


27.41 


2, 


.52 


16073, 


.7601 


29.82 


2 


.49 



Table 5 

Radial Velocities for HD 71881 



HJD RV Uncertainty 

-2,440,000 (m s^i) (m s''^) 



10807, 


.1721 


112.06 


2, 


.56 


10837, 


,9895 


106.74 


2, 


.37 


10862, 


.7614 


109.55 


2, 


.46 


10862, 


.8804 


108.43 


2, 


.37 


11171, 


.0933 


97.73 


2 


.35 


11227, 


.9372 


98.27 


2, 


.32 


11552, 


.0265 


79.55 


2, 


.42 


11900, 


.1047 


76.05 


2 


.33 


12003, 


.9078 


69.54 


2, 


.47 


12007, 


.8729 


66.09 


2, 


.39 


12062, 


.7530 


69.57 


2 


.47 


12064, 


.7717 


71.06 


2, 


.55 


12236, 


.0510 


56.62 


2 


.54 


12573, 


.1159 


51.72 


2 


.40 


12712, 


.8234 


48.80 


2, 


.50 


12988, 


.0181 


45.16 


2, 


.43 


13303, 


.1333 


50.16 


3 


.91 


13425, 


.8843 


25.04 


4, 


.67 


15556, 


.1133 


-8.53 


4, 


.41 


15871, 


.0606 


-34.85 


4, 


.41 


KiOLS, 


,8872 


-40.97 


4 
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